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Abstract. We will present a summary of the observations of blazars by the 
Energetic Gamma Ray Experiment Telescope (EGRET) on the Compton 
Gamma Ray Observatory (CGRO). EGRET has detected high energy 7- 
ray emission at energies greater than 100 MeV from more that 50 blazars. 
These sources show inferred isotropic luminosities as large as 3 x 

1Q 49 

eras 

s _1 . One of the most remarkable characteristics of the EGRET observations 
is that the 7-ray luminosity often dominates the bolometric power of the 
blazar. A few of the blazars are seen to exhibit variability on very short 
time-scales of one day or less. The combination of high luminosities and 
time variations seen in the 7-ray data indicate that 7-rays are an important 
component of the relativistic jet thought to characterize blazars. Currently 
most models for blazars involve a beaming scenario. In leptonic models, 
where electrons are the primary accelerated particles, 7-ray emission is 
believed to be due to inverse Compton scattering of low energy photons, 
although opinions differ as to the source of the soft photons. Hardronic 
models involve secondary production or photomeson production followed 
by pair cascades, and predict associated neutrino production. 



1. Introduction 

One of the most striking accomplishments of the Energetic Gamma Ray 
Experiment Telescope (EGRET) instrument on the Compton Gamma-Ray 
Observatory (CGRO) is the detection of high-energy 7-rays from active 
galaxies whose emission at most wavebands is dominated by non-thermal 

* Invited review paper to appear in Observational Evidence for Black Holes in the 
Universe, ed. S. K. Chakrabarti, 1999 (Dordrecht: Kluwer) pg 215-230. 
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processes. These objects, called "blazars," are highly variable at most fre- 
quencies and are bright radio sources. Prior to the launch of CGRO, 3C 
273, discovered by COS-B (Swanenburg et al. 1978), was the only known 
extragalactic source of 7-rays. Since then, EGRET has detected more than 
50 blazars in high energy (> 100 MeV) 7-rays (Mukherjee et al. 1997; 
Thompson et al. 1995; 1996). 

The blazars detected by EGRET all share the common characteristic 
that they are radio-loud, flat-spectrum radio sources, with radio spectral 
indices a r > —0.6 (von Montigny et al. 1995). Several of these blazars are 
known to demonstrate superluminal motion of components resolved with 
VLBI (3C 279, 3C 273, 3C 454.3, PKS 0528+134, for example). The blazar 
class of active galactic nuclei (AGN) includes BL Lac objects, highly polar- 
ized quasars (HPQ), or optically violent variable (OVV) quasars and arc 
characterized by one or more of the properties of this source class, namely, 
a non-thermal continuum spectrum, a flat radio spectrum, strong variabil- 
ity and optical polarization. For many of the EGRET-detected blazars, the 
7-ray energy flux is dominant over the flux in lower energy bands. The 
redshifts of these sources range from 0.03 to 2.28 and the average photon 
spectral index, assuming a simple power law fit to the spectrum, is ~ 2.2. 
Many of the blazars exhibit variability in their 7-ray flux on timescales 
of several days to months. In addition, blazars exhibit strong and rapid 
variability in both optical and radio wavelengths. 

Of the 51 blazars reviewed here, 14 are BL Lac objects, and the rest 
are flat spectrum radio quasars (FSRQs). BL Lac objects generally have 
stronger polarization and weaker optical lines. In fact, some BL Lac ob- 
jects have no redshift determination because they have no identified lines 
above their optical continuum. FSRQs are generally more distant and more 
luminous compared to the BL Lac objects. 

This review summarizes the present knowledge on 7-ray observations of 
blazars by EGRET. A brief description of the EGRET instrument and data 
analysis techniques, and the list of blazars detected by EGRET is given in 
§2. Temporal variations and 7-ray luminosity of blazars are discussed in 
§§3 & 4. Section 5 describes the spectral energy distribution of blazars and 
summarizes the various models that have been proposed to explain the 
7-ray emission in blazars. 

2. EGRET observations and analysis 

2.1. THE EGRET INSTRUMENT 

EGRET is a 7-ray telescope that is sensitive in the energy range ~ 30 MeV 
to 30 GeV. It has the standard components of a high-energy 7-ray instru- 
ment: an anticoincidence dome to discriminate against charged particles, 
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a spark chamber particle track detector with interspersed high-Z material 
to convert the 7-rays into electron-positron pairs, a triggering telescope to 
detect the presence of the pair with the correct direction of motion, and 
an energy measurement system, which in the case of EGRET is a Nal(Tl) 
crystal. EGRET has an effective area of 1500 cm 2 in the energy range 0.2 
GeV to 1 GeV, decreasing to about one-half the on-axis value at 18° off-axis 
and to one-sixth at 30° . The instrument is described in details by Hughes et 
al. (1980) and Kanbach et al. (1988, 1989) and the prefhght and postflight 
calibrations are given by Thompson et al. (1993) and Esposito et al. (1998), 
respectively. 

Although EGRET records individual photons in the energy range 30 
MeV to about 30 GeV, there are several instrumental characteristics that 
limit the energy range for which time variation investigations of blazars 
are viable. At the low end of the energy range, below ~ 70 MeV, there are 
systematic uncertainties that make the spectral information marginally use- 
ful. In addition, the deteriorating point spread function (PSF) and energy 
resolution at low energies, make analysis more difficult. At high energies, 
although the systematic uncertainties are reduced, and the PSF and energy 
resolution are more reasonable, because of the steeply falling spectra, few 
photons are detected above 5 GeV. 

The angular resolution of EGRET is energy dependent, varying from 
about 8° at 60 MeV to 0.4° above 3 GeV (68% containment). The positions 
of sources are detected with varying accuracy: better than 0.1° for the very 
bright sources, or at least 0.5° for sources just above the detection threshold. 

The threshold sensitivity of EGRET (> 100 MeV) for a single observa- 
tion is ~ 3 x 10~ 7 photons cm~ 2 s _1 , and is only about a factor of 50-100 
below the maximum blazar flux ever observed. The dynamic range for most 
observations of blazar variations is, therefore, fairly small. 

2.2. EGRET DATA ANALYSIS 

The blazars described here were typically observed by EGRET for a period 
of 1 to 2 weeks; however, several of them were observed for 3 to 5.5 weeks. 
Following the standard EGRET processing of individual 7-ray events, sum- 
mary event files were produced with 7-ray arrival times, directions and 
energies. For the observations reported here, photons coming from direc- 
tions greater than 30° from the center of the field of view (FOV) were not 
used, in order to restrict the analysis to photons with the best energy and 
position determinations. In addition, exposure history files were produced 
containing information on the instrument's mode of operation and point- 
ing. These maps were used to generate skymaps of counts and intensity for 
the entire field of view for each observation, using a grid of 0.5° x 0.5°. The 
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intensity maps were derived simply by dividing the counts by the exposure. 
The EGRET data processing techniques are described further by Bertsch 
et al. (1989). 

The number of source photons, distributed according to the instrument 
PSF in excess of the diffuse background, was optimized. An E~ 2 photon 
spectrum was initially assumed for the source search. The background dif- 
fuse radiation was taken to be a combination of a Galactic component 
caused by cosmic ray interactions in atomic and molecular hydrogen gas 
(Hunter et al. 1997), as well as an almost uniformly distributed component 
that is believed to be of extragalactic origin (Sreekumar et al. 1998). 

The data were analyzed using the method of maximum likelihood as de- 
scribed by Mattox et al. (1996) and Esposito et al. (1998). The likelihood 
value, L, for a model of the number of 7-rays in each pixel of a region of the 
map is given by the product of the probability that the measured counts 
are consistent with the model counts assuming a Poisson distribution. The 
probability of one model with likelihood, L\, better representing the data 
than another model with likelihood, L2, is determined from twice the dif- 
ference of the logarithms of the likelihoods, 2(lnL2 — InLi). This difference, 
referred to as the test statistic TS, is distributed like % 2 with the number 
of degrees of freedom being the difference in the number of free parameters 
in the two models. The flux of the point source and the flux of the diffuse 
background emission in the model are adjusted to maximize the likelihood. 
The significance of a source detection in sigma is given approximately by 
the square root of TS. 



2.3. EGRET OBSERVATIONS 



The 51 blazars listed in Table 1 were all detected by EGRET above 100 
MeV during the period of EGRET observations from 1991 April to 1995 
September (Phases 1 through 4 of CGRO) (Mukherjee et al. 1997). Some of 
these blazar associations are not certain; Mattox et al. (1997a) find only 42 
identifications to have high confidence. Conversely, some of the unidentified 
high-latitude EGRET sources are likely to be blazars. In addition to the 
42 considered strongest, Mattox et al. (1997a) note 16 possible associations 
with bright flat-spectrum, blazar-like radio sources. Typically, each blazar 
listed in Table 1 was seen in several different viewing periods (VPs). The 
maximum and minimum fluxes observed for each blazar is indicated in 
Table 1. A more complete list of blazar detections by EGRET may be 
found in the third EGRET catalog (Hartman et al. 1998). 



High Energy Gamma-Ray Emission from Blazars: 



5 



TABLE 1 

EGRET-DETECTED BLAZARS 
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0.84 
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a Redshift. 

b Units: 10 -8 ph cm -2 s -1 . Upper limits are 2 sigma. 
c Variability index (see text for definition). 

d Both the maximum flux and the variability index would increase if the recent flare in BL Lac were 
taken into account (Bloom et al. 1997). 
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Figure 1. Flux history of four blazars (PKS 1633+382, PKS 1406-076, Mrk 421, and 
PKS 0528+134) from 1991 April to 1995 September, as detected by EGRET. 2a upper 
limits are indicated by downward arrows. 



3. Time variability 

The fluxes of the blazars detected by EGRET have been found to be vari- 
able on time scales of a year or more down to well under a day. Long term 
variations of blazars have been addressed earlier by several authors (eg. 
von Montigny et al. 1995; Hartman et al. 1996a; Mukherjee et al. 1997). In 
some cases, many of the detected blazars have exhibited flux variations up 
to a factor of about 30 between different observations. Figure 1 shows the 
flux history of four EGRET-detected blazars. The horizontal bars on the 
individual data points denote the extent of the VP for that observation. 
Fluxes have been plotted for all detections greater than 2a. For detections 
below 2<7, upper limits at the 95% confidence level are shown. A systematic 
uncertainty of 6% was added in quadrature with the statistical uncertainty 
for each flux value, consistent with the analysis of McLaughlin et al. (1996) 
on EGRET source variability. 

In order to quantify the flux variability of the blazars in Table 1 , Mukher- 
jee et al. (1997) calculated the variability index, V = logQ, as defined by 
McLaughlin et al. (1996), where Q = 1 — P x (x 2 ,v)- Here P x (x 2 ,v) is the 
probability of observing x 2 ° r something larger from a x 2 distribution with 
v degrees of freedom. The flux versus time data were fit to a constant flux 
and the reduced xli f° r v degrees of freedom, was calculated using the least 
square fit method. For a nonvariable source, a constant flux is expected to 
fit the data well, and the mean value of the x 2 distribution is expected to 
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Figure 2. Distributions of the variability indices for BL Lac objects and FSRQs. 



be equal to the number of degrees of freedom in the data. The quantity V 
is used to judge the strength of the evidence for flux variability. Following 
the classification of McLaughlin et al. , V < 0.5 was taken to indicate non- 
variability, V > 1 to indicate variability, and 0.5 < V < 1 as uncertain. 
Table 1 lists the value of V for each source. 

Of the 51 blazars reviewed here, 35 are found to be variable (V > 1), 9 
are non- variable (V < 0.5), and 7 fall in the range of uncertain variability. It 
should be noted that, although the criterion used here to gauge variability 
of a blazar is somewhat arbitrary, it does provide a way to compare the 
numbers obtained. Also, as McLaughlin et al. (1996) have noted, changing 
these criteria by 20% yields similar results. If the FSRQs and the BL Lac 
objects are considered separately, it is found that 76% of the FSRQs in the 
sample are variable, while 16% are non-variable. Similarly, for the BL Lac 
objects, 50% are variable, while 21% are definitely non-variable. It should 
be noted that the low intrinsic luminosity of BL Lac objects could bias 
observations (see discussions in §5.2). Figure 2 shows the distribution of 
the variability indices for the FSRQs and the BL Lac objects. The BL Lac 
objects in the data set are found to be less variable on the average than the 
FSRQs. Recent observations of flares in BL Lac objects, however, modify 
some of these conclusions. For example, BL Lac was detected during a 7- 
ray outburst with an average flux of (171 ± 42) x 10~ 8 photons cm~ 2 s -1 
in July 1997 (Bloom et al. 1997). BL Lac would have a high value of V in 
Table 1, if this information were taken into account. 

Figure 3 shows a plot of the variability index as a function of the 
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Figure 3. Variability index plotted as a function of the weighted average flux for BL 
Lac objects and FSRQs. The BL Lac objects are denoted by diamonds and the FSRQs 
by plus signs in the figure. 

weighted average flux for the blazars in Table 1. Note that the sources 
that have the highest average fluxes all have high variability indices. Only 
2 out of 18 blazars with average flux less than 1 x 10~ 7 photons cm -2 s _1 
have a variability index greater than 2.5. In fact, there are no non-variable 
blazars in the sample that have high flux. 

The study of short term variability in blazars is always limited by the 
small numbers of photons detected in the short time intervals at 7-ray ener- 
gies. The shortest time-scale variations detected for blazars with EGRET 
are for PKS 1622-297 (Mattox et al. 1997b) and 3C 279 (Wehrle et al. 
1997). For both these objects the flux was found to increase by a factor 
of two or more in less than 8 hours. Other objects that have shown flux 
variations over the period of a few days are 3C 279 (Kniffen et al. 1993), 3C 
454.3 (Hartman et al. 1993), 4C 38.41 (Mattox et al. 1993), PKS 1406-076 
(Wagner et al. 1995), and PKS 0528+134 (Hunter et al. 1993; Mukherjee 
et al. 1996). The short time-scale of 7-ray flux variability (e.g. in 3C 279 
or PKS 1622-297) when combined with the large inferred 7-ray luminosi- 
ties, implies that the blazar emission region is very compact. Gamma-ray 
tests for beaming from variability and flux measurements using the Elliott- 
Shapiro relation and 7-ray transparency arguments are summarized in a 
recent review on 7-ray blazars by Hartman et al. (1997). A factor-of-two 
flux variation on an observed time-scale <5t bs limits the size r of a station- 
ary isotropically emitting region to be roughly r < c5t ^, s /(l + z) by simple 
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light-travel time arguments. Under the assumptions of isotropic radiation 
and Eddington-limited accretion, the implied minimum black hole masses 
of blazars are > 8 x 10 11 M for PKS 1622-297 (Mattox et al. 1997b) 
from EGRET observations and > 7.5 x 10 s M Q for PKS 0528+134 from 
COMPTEL observations (Collmar et al. 1997). 

4. Luminosity 

The 7-ray luminosity can be estimated by considering the relationship be- 
tween the observed differential energy flux Sq(Eq), where the subscript "0" 
denotes the observed or present value, and Q e the power emitted in dE, 
where E = Eq(1 + z) in the Friedman universe. 

Q e [E ] = 47tS (£o)(1 + z^QDl^z, q ) (1) 

where 

°L = 17—2^ " 90 + qoz + (q - l)(2q z + l) 1 / 2 ] = ^-g(z, q ). (2) 

Hq is the Hubble parameter, q$ is the deceleration parameter, b is the 
spectral index, z is the redshift, and O is the beaming factor. Hq is chosen 
to be 70 and qo to be 0.5, although the results obtained here are not highly 
sensitive to these choices. The beaming factor is taken to be 1. The spectral 
index is obtained using the analysis described in §4. The luminosity as a 
function of the redshift is determined using equation (1) and is plotted 
in Figure 4 for the blazars detected by EGRET. The typical detection 
threshold for EGRET as a function of z, for relatively good conditions, 
is also shown in the figure. The actual threshold varies somewhat with 
exposure and region of the sky, and the average threshold is a little higher 
than the curve shown, but the shape is the same. The BL Lac objects are 
indicated in the figure with dark diamonds, and one sees clearly that they 
are predominantly closer and lower in luminosity. 

Recently, Chiang & Mukherjee (1998) have calculated the evolution and 
luminosity function of the EGRET blazars, and have estimated the con- 
tribution of this source class to the diffuse extragalactic gamma-ray back- 
ground. They find that the evolution is consistent with pure luminosity evo- 
lution. According to their estimates, only 25% of the diffuse extragalactic 
emission measured by SAS-2 and EGRET can be attributed to unresolved 
7-ray blazars, contrary to some of the other estimates (eg. Stecker & Sala- 
mon 1996). Below 10 MeV, the average blazar spectrum suggests that only 
about 50% of the measured 7-ray emission could arise from blazars (Sreeku- 
mar, Stecker & Kappadath 1997). This leads to the exciting possibility that 
other sources of diffuse extragalactic 7-ray emission must exist. 
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Figure 4- Luminosity vs redshift for blazars detected by EGRET. The BL Lac objects 
are indicated with filled symbols. The typical detection threshold for EGRET is shown 
as a solid curve. 

5. Spectra 

5.1. SPECTRA IN THE EGRET ENERGY RANGE 

EGRET spectra of blazars typically covers at least two decades in energy 
(from 30 MeV to 10 GeV) and are well described by a simple power-law 
model of the form F(E) = k(E/Eo)~ a photons cm~ 2 s _1 McV -1 , where 
the photon spectral index, a, and the coefficient, k, are the free parameters. 
The energy normalization factor, Eq, is chosen so that the statistical errors 
in the power law index and the overall normalization are uncorrelated. 

The average blazar spectrum has a spectral index of about —2.15. Figure 
5 shows the photon spectral index of the blazars plotted as a function of 
the redshift. There are marginal indications that suggest that the BL Lac 
objects have slightly harder spectrum in the EGRET energy range than the 
FSRQs. Mukherjee et al. (1997) find that the average spectral index of the 
BL Lac objects is about —2.03, compared to about —2.20 for the FSRQs. 
For some individual blazars there has been noted a trend for the spectrum 
to harden during a flare state (eg. in blazars 1222+216, 1633+382, and 
0528+134; Sreekumar et al. 1996; Mukherjee et al. 1996). A spectral study 
of blazars as a class has been performed by Miicke et al. (1996) and Pohl 
et al. (1997). 
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Figure 5. Photon spectral index as a function of the redshift for blazars detected by 
EGRET. 

5.2. SPECTRAL ENERGY DISTRIBUTIONS AND GAMMA-RAY MODELS 

The processes by which 7-rays are produced in blazars can be best under- 
stood by the study of the correlated multiwaveband observations of blazars 
extending from radio to 7-ray wavebands. One of the most significant find- 
ings of EGRET is that, in the radio to 7-ray multiwavelength spectra of 
blazars, the power in the 7-ray range equals or exceeds the power in the 
infrared-optical band. Any model of high-energy 7-ray emission in blazars 
needs to explain this basic observational fact. The high 7-ray luminosity of 
the blazars suggests that the emission is likely to be beamed and, therefore, 
Doppler-boosted into the line of sight. This is in agreement with the strong 
association of EGRET blazars with radio-loud flat-spectrum radio sources, 
with many of them showing super luminal motion in their jets. This infor- 
mation has helped to favor jet models of emission over models in which the 
7-ray production is directly associated with accretion onto a massive black 
hole (e.g. Becker & Kafatos 1993). 

The jet models explain the radio to UV continuum from blazars as 
synchrotron radiation from high energy electrons in a relativistically out- 
flowing jet which has been ejected from an accreting supermassive black 
hole (Blandford & Konigl 1979). The emission in the MeV-GeV range is 
believed to be due to the inverse Compton scattering of low-energy photons 
by the same relativistic electrons in the jet. However, two main issues re- 
main questionable: the source of the soft photons that are inverse Compton 
scattered, and the structure of the inner jet, which cannot be imaged di- 
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rectly. The soft photons can originate as synchrotron emission either from 
within the jet (the synchrotron-self-Compton or SSC process: Maraschi, 
Ghisellini, & Celotti 1992; Bloom & Marscher 1996), or from a nearby 
accretion disk, or they can be disk radiation reprocessed in broad-emission- 
line clouds (the external radiation Compton process or the ERC process: 
Dermer & Schlickeiser 1994; Sikora, Begelman, & Rees 1994; Blandford & 
Levison 1995; Ghisellini & Madau 1996). In contrast to these leptonic jet 
models, the proton-initiated cascade (PIC) model (Mannheim & Biermann 
1989, 1992) predicts that the high-energy emission comes from knots in jets 
as a consequence of diffusive shock acceleration of protons to energies so 
high that the threshold of secondary particle production is exceeded. 

Figure 6 shows the simultaneous spectral energy distribution of 3C 279 
during January-February 1996, when the source was detected at its highest 
state ever (Wehrle et al. 1997). The figure shows the relative amounts of 
energy detected in equal logarithmic frequency bands. The power output 
in 7-rays dominates the bolometric luminosity of the sources, as mentioned 
in §1. Wehrle et al. (1997) note that the 7-rays vary by more than the 
square of the observed IR-optical flux change, a fact that could be hard to 
explain by some specific blazar emission models. Although the data do not 
rule out SSC models, Wehrle et al. point out that the data are most likely 
explained by the "mirror" model of (Ghisellini & Madau 1996). In this 
model the flaring region in the jet photo- ionizes nearby broad-emission- line 
clouds, which in turn provide low energy external seed photons that are 
inverse Compton-scattered to high energy 7-rays. 

Recently, a model combining the ERC and SSC scenarios has been used 
to fit the simultaneous COMPTEL and EGRET spectra of PKS 0528+134 
by Bottcher & Collmar (1998). Figure 7 shows their fit to the gamma-ray 
spectrum during the high state of the source during March 1993. In their 
model Bottcher &: Collmar assume a spherical blob filled with ultrarelativis- 
tic pair plasma which is moving out along an existing jet structure perpen- 
dicular to an accretion disk around a black hole of mass 5 x 10 10 M . They 
argue that the observed spectral break between COMPTEL and EGRET 
energy ranges can plausibly be explained by a variation of the Doppler 
beaming factor in the framework of a relativistic jet model for AGNs. 

The EGRET results have demonstrated that in order to model the spec- 
tra of blazars it is very important to get a truly simultaneous coverage across 
the entire electromagnetic spectrum before, during, and after a flare in the 
high-energy 7-ray emission. The limited data that we have on most of the 
blazars prevents us from being able to distinguish between the different 
theoretical models, on the basis of the spectra alone. For example, both the 
SSC and ERC models have been shown to reproduce the multiwavelength 
spectrum of 3C 279 rather well (Hartman et al. 1996b; Maraschi, Ghisellini 
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Figure 6. Radio to 7-ray energy distribution of 3C 279 in low (open circles) and flaring 
state (filled circles) in f 996 January- February (Wehrle et al. f 997) . 



& Celotti 1992; Ghisellini & Maraschi 1996). The SSC model was similarly 
found to fit the multiwavelength spectrum of PKS 0528+134 during the 
March 1993 flare reasonably well (Mukhcrjce et al. 1996). The low-state 
data of PKS 0528+134 (Aug 1994) was fit well with the ERC model, as 
demonstrated by Sambruna et al. (1997). The SSC, ERC, and PIC models 
have all been shown to fit the multiwavelength spectrum of 3C 273 well 
(von Montigny 1997). 

The differences between the 7-ray variability properties of BL Lac ob- 
jects and FSRQs can be explained in light of the model of Ghisselini and 
Madau (GM) (1996). In their model, soft photons from the jet are repro- 
cessed by broad line region (BLR) clouds. Subsequently, these soft photons 
are emitted back into the jet where they scatter off of electrons in relativis- 
tically moving "blobs" to create high-energy 7-rays. Since BL Lac objects 
generally have very weak emission lines, it may be that they have much 
less BLR gas available for reprocessing soft photons. If there is an initial 
outburst of soft photons created via the synchrotron process in a jet "blob," 
then BL Lac objects can still create 7-rays via the SSC process, though per- 
haps with lower amplitude than 7-rays created via the GM model. (This 
effect is somewhat dependent on adjustable model parameters.) In this sce- 
nario, BL Lac objects may undergo several SSC outbursts which fail to 
reach the EGRET detection threshold, thus giving the appearance that BL 
Lac objects as a class experience less dramatic and less frequent 7-ray flares. 
The general properties of the low frequency outbursts of BL Lac objects, 
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Figure 7. Fit to the 7-ray (COMPTEL and EGRET) spectrum of PKS 0528+134 in 
its high state in VP 213. See Bottcher & Collmar (1998) for fit parameters. 

however, would be very similar to that of the FSRQs. 

In order to achieve a better understanding of the emission mechanisms of 
7-rays from blazars a study of the correlated short time scale (~ 1 — 3 days) 
7-ray variations with those at other frequency bands is needed. Since the 
predictions of time delays between the flux changes at various frequencies 
are different for the individual models for both the seed photons and the 
nature of the inner jet, this method could provide a means to discriminate 
between the different models. The differences in the model predictions are 
discussed in more detail by Marscher et al. (1995). Gamma-ray variability 
in the different models may have different impacts on the spectral behavior 
during the build-up and decline of an outburst. Studying the short-time- 
scale behavior and looking for spectral changes while following a complete 
outburst may be the key to pin down the basic emission mechanisms. 

6. Summary 

In conclusion, the EGRET results have shown the importance of the 7- 
ray window on blazars. The high luminosities and strong time variability 
observed have pushed theoretical models to emphasize relativistic jets of 
particles seen at small angles to the line of sight. The EGRET observations 
have established that the 7-ray window is critical for understanding the 
properties of blazars. Future observations with CGRO and successor 7- 
ray observatories like INTEGRAL and GLAST should play a key role in 
resolving the physics of these powerful sources. 
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